Interaction effects on random Dirac fermions 
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We study a Dirac fermion model with three kinds of disorder as well as a marginal interaction 
which forms the critical line of c = 1 conformal field theory. Computing scaling equations by the 
use of a perturbative renormalization group method, we investigate how such an interaction affects 
the universality classes of disordered systems with non-interacting fermions. We show that some 
specific fixed points are stable against an interaction, whereas others are unstable and flow to new 
random critical points with a finite interaction. 
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I. INTRODUCTION 



Anderson localization has attracted renewed interest in 
recent years, especially due to the discovery by Altland 
and Zirnbaueru of novel universality classes for disordered 
systems. In addition to the well-known Wigner-Dyson 
classification of the random matrices, seven new classes 
have been established. They can be found in quasipar- 
ticle systems described by Bogoliubov-de Gennes Hamil- 
tonians, or in systems with sublattice (chiral) symmetry, 
studied for the first time by Gade and Wegner.0 As is 
well-known, symmetry properties play a crucial role in 
the localization problems. The key symmetry of the new 
classes is a particle-hole symmetry, which makes the spec- 
trum symmetric with respect to the zero-energy. 

The new universality classes have turned out to be 
quite useful in studying problems of quasinarticle lo- 
calization in dirty d-wave superconductorsJ3i3 Random- 
bond models on bipartite lattices such as the random 
flux modelu-the random hopping model with it flux per 
plaquettejj113 graphite sheets,E3 etc have also been ex- 
amined successfully from the point of view of the chi- 
ral classes. What is interesting is that many of these 
models can be described, near their critical points, by 
Dirac fermions with appropriate symmetries. Surpris- 
ingly, Dirac fermion models, with a slight extension 
to multi-species cases, exhaust all possible universality 
classes found so far, and detailed classification has been 
given by Bernard and LeClairO 

On the other hand, a disordered Dirac fermion model 
has been proposed as an .effective model for the integer 
quantum Hall transition l 14 i Actually, the Chalker and 
Coddington network modeta, which is believed to de- 
scribe the transition correctly, yields disordered Dirac 
Hamiltonian in a continuum limit $3 Some specific fixed 
points of the mO|dfjL,e^f;. the fixed line due to a random 
vector potj^tial,llJO't3 or the fixed point due to a ran- 
dom masstJcj have been analyzed in a weak-coupling 
regime. 

Although the disordered Dirac fermions are one of the 
simplest models, many problems still remain to be ex- 
plored. One of them is to clarify strong coupling regime of 



the model. It is quite necessary to find out a generic fixed 
point of the Dirac fermions with generic disorder, because 
it should be a field theory describing the integer quantum 
Hall transition. Moreover, even the well-known fixed line 
formed by a random vector potential should be in a dif- 
ferent phase, in a strong coupling regime, from the|-aMj 
expected by conventional weak-coupling approaches. r 3 lr 4 l 
However, recent developments enable us to discuss strong 
coupling aspects by the use of thr ijruormalization group 
method or a variational method.E3iHj 

Another interesting problem is to clarify effects of in- 
teraction on the localization properties. Studies on the 
random bond Ashkin-Tcllcr._model have been reported 
by Dotsenko and DotsenkorH Describing the model by 
Majorana fermions and computing one loop renormal- 
ization group (RG) equations, they have found that the 
Harris criteria does not hold any longer in this model. 
Their question is whether this feature is common in two- 
dimensional systems. More generally, an interplay be- 
tween interactions and randomness has been extensively 
studiedB3 

The notion of the new universality classes, mentioned- 
above, for the quasiparticle systems as well as random 
hopping models on bipartite lattices are only for non- 
interacting systems, and if models include interactions 
and randomness simultaneously, interactions in general 
break the symmetry which the models without interac- 
tions should have. Therefore, it is quite interesting to 
study how interactions affect the symmetry properties of 
the newly discovered universality classes. 

In this paper, we study a Dirac fermion model with 
three kinds of disorder as well as a marginal interaction 
which forms the critical line of the c = 1 conformal field 
theory, in order to clarify whether the universality classes 
for non-interacting theories are stable against such an in- 
teraction. Conversely, we might say that we study how 
c = 1 theory are changed by randomness. We calculate 
one-loop RG equations by the use of the replica trick. 
It turns out that there is an interesting interplay be- 
tween the marginal interaction and randomness. Namely, 
the universality classes for non-interacting fermions are 
stable again this interaction in some cases, but in other 
cases, the interaction survives and show the critical line 



1 



like the c = 1 theory. 

This paper is organized as follows. In the next section, 
we define the model and derive the scaling equations. 
From Sec. Ill to Sec. VI, we study four kinds of specific 
fixed points in detail. In Sec. VII, we summarize the 
results and briefly discuss a generic fixed point. 



II. THE MODEL AND SCALING EQUATIONS 

The model we will study in this paper is described by 
the following Dirac Hamiltonian in two dimensions with 
quenched disorder as well as an interaction, 



= Iphlp + Hint, 



(2.1) 



where M(x), A fl (x), and V(x) are a random mass, a 
random vector potential, and a random scalar potential, 
respectively, and ip and if) fields are two-component Dirac 
fields, whose chiral components are defined by 

With respect to these notations the interaction Tti n t is 
defined here as 



'Hint — 2gJn,JL, 
where R- and L-currents are defined by 



(2.3) 



(2.4) 



Without disorder, this interaction is exactly marginal, 
and the model is on the critical line with continuously 
varying critical exponents of the c — 1 conformal field 
theory. 



The model (2J) can describe critical behavior of var- 
ious kinds of lattice models. One of well-known exam- 
ples is the Ashkin- Teller model or equivalently the iBafc 
ter model. Without randomness, it has been shownC3Lll 
that the models form a critical line, and in the scaling 
limit, they can be described by two kinds of Majorana 
fermions or eq uivalently one Dirac fermion with the in- 
teraction (2.3). The case g = corresponds to the decou- 
pling point of two Ising spins, giving just two independent 
Ising models. 

If random bond interactions are introduced to this 
model, the Dir ac fermion acquires a ranclpm mass term 
as in Eq. (2J). Dotsenko and Dotsenkdij have firstly 
studied the Ising model (the g — case) and shown that 
the random mass is marginally irrelevant, and it gives 
rise to the famous In | lnr| behavior of the specific heat, 
where r is the reduced temperature. In the case.of the fi- 
nite interaction, the same authors have showncH that the 
positive interaction leads to the Ising fixed point g = 0, 
whereas the negative interaction leads to l/(ln | lnr|) be- 
havior of the specific heat. This means that interaction 



and randomness couple together, yielding new critical be- 
havior. 

Other lattice models are two-dimensional quantum sys- 
tems of (nonrelativistic) fermion on bipartite lattices such 
as a square lattice with 7r-flux per plaquette and a honey- 
comb latt ice. In such quantum systems, the interaction 
term (2.3) does not arise from interactions of the origi- 
nal lattice models, since one has to treat, in that case, 
a field theory of the Dirac fermion in three dimensions. 
One possibility is the coexistence of annealed randomness 
and quenche d r andomness: The former can serve as the 
interaction (2.3). Actually, Belitz et. al. have stressed 
the importance of annealed randomness for systems with 
quenched randomness and proposed a new mechanism 
for a metal-insulator transition.ES 

The partition function is defined by 



Z 



VipVip exp 



d 2 3 



H 



(2.5) 



where z = E + ie. In order to compare the aspects of sev- 
eral fixed points, we will devote ourselves to the calcula- 
tion of a two-point function and resultant density of state 
(DOS). For this purpose, we int roduced a source term 
zipip i n th e action functional ( |2.5| ). Without the interac- 
tion (2^3), the two-point f unct ion is {ipip) = (h — z)^ 1 , 
where h is defined in Eq. fl2.ip . In the present case, the 
interaction Timt yields a self-energy £ and the two-point 
function can be expressed as (iptp) — (h — £ — z)~ . 

Now we assume that the three kinds of randomness 
obey the Gaussian distribution with zero mean 



A^{x)A v {y) = gAS^S{x - y), 
M(x)M(y)=g M 6(x-y), 
V(x)V{y) = g v S(x - y). 



(2.6) 



Then ensemble-average over disorder by the use of the 
replica trick leads to the following Hamiltonian, 



H = ip a ijudu-p a ~ y {^alnipaY - (ip a <J3^ a y 

(V'qV'q) 2 + ZgJRaJLa ~ Z1p a 1p a . (2.7) 



s 2 gn 



where a = 1,2, 3.. .n denote the replicas, and Jr(l)ci = 

a • 

One can compute the 1-loop beta functions using the 
operator product expansion. The unperturbed two-point 
functions read, 



(ipR<x(z)ip Ra i(w)) = 
(tpL a (z)ip La ,(w)) = 



3aa' 1 

2n z — w' 



2tt z — w 

and using these, we have the scaling equations 



(2.8) 



2 



d 9A 1 / 2 2\ 

a ^- = -[Ag A + {l-n)g_+2g]g + , 
dl 7r 



dl 2tt LV 



ng_ 



4gg-\ 



1 

dl 7T 

1 dz _ ! 

2 oH 



— [2 SA + (1 - n) 5+ - n<?_ + 2g] , (2.9) 

Z7T 



where g± = gyigM, and / = InL with the system length 
L. Taking the replica limit n — > 0, we reach the follow- 
ing RG equations which we will study in the following 
sections, 



dl 



~9% 



= 9l 



99- 



dg+ ,~ ~ ~\ ~ 

= {9 a + 9- +9)9+, 

dg- 
dl 
dg 

1 dz „ 
z dl 



(2.10) 



where we have denoted the coupling constants as cja — 
^9a/-k, g± = g±/ir, and g = 2g/w. When g = 0, these, 
equations are the same as those derived by Bernard.Ea 
The anomalous dimension 7 of ipip is 



1 1 9a .9+ , g 

7 = 1— — 

1 4 2 2 



(2.11) 



The last equation in ( 2.1C| ) tells us that the energy 
E = Rez grows according to 



A r l 

In- = 2(/-I )-y j 



(2.12) 



where E denotes the bare energy at the length scale 
Iq = lnL while A is the renormalized energy at the 
scale I = InL and serves as a cut-off. Since the scal- 



ing dimension of the field iptp is given by (2.11), the 
ensemble-averaged DOS, which is the imaginary part of 
lim e ^o tr^i/j), should obey the scaling law, 



P(L/L ) 



exp 



dl-f 



(2.13) 



Ludwig et.|_al. have studied the same model without 
interaction.tj They examined several specific fixed points 
of the model, although in the generic case with all pos- 
sible randomness, scaling equations lead to a strong cou- 
pling regime, which is not accessible by a perturbative 
weak-coupling approach. In the following sections, we 
will examine how several fi xed points known so far are 
affected by the interaction (2.3). 



III. THE RANDOM VECTOR POTENTIAL 

In this section, we study the model with the random 
vector potential and the interaction only. First of all, it 
should be noted that without the interaction, the Hamil- 
tonian has chiral U(n)xU(ji) symmetry, and therefore, 
belongs to the class AIII.El The interaction, however, 
breaks it up to n copies of U(l)xU(l) symmetry. 

It is known that the random vector potential is exactly 
marginal. On the other hand, the interaction is also ex- 
actly marginal. It turns out that even if they exist at the 
same time, they remain to be exactly marginal. Actually, 
by setting g + = g_ = 0, the scaling equations reduce to 



d£A _ a jj_ = n 
dl dl 



(3.1) 



Therefore, the anomalous dimension is a constant given 
by 



= 1 _ - 91 _ 1 _ 9 AO + go 

1 ~ 4 2 7T 



(3.2) 



where gAo = ^gAo/'K and go = .90/ 71 " ar< 3 initial coupling 
constants at I = Iq. The energy flow is 



, A , . gAo + go 

and hence the density of state is given by 

P{E) ~ (a) • 



(3.3) 



(3.4) 



When go =^S) this formula is just the same as that by Lud- 
wig et. alli-l The interaction just modifies the exponent. 
However, from the point of view of symmetry properties, 
the model with finite interaction does not belong to AIII 
any longer. 



IV. THE RANDOM MASS 

It is well-known that the random-bond Ising model can 
be described, near the transition temperature, by the 
present Dira c model with the random mass term only 
in Eq. (2.1).EfE3 This model is invariant under 0(2n) 
rotation with respect to Maiorana fermions, and there- 
fore belongs to the class The Ashkin- Teller model, 
which is composed of two Ising models coupled with 4- 
spin interaction, 



H = - >^ [Jij(o-ii<rij + cr 2 iCT2j) + J4(J ll ai j a2ia 2 j} , (4.1) 



(id) 



can be described by the Dirac model with the interac- 
tion (2.3), and with thej-^jandom mass as well, if bond- 
randomness is included.EII As in the previous case, the 



3 



interaction breaks the symmetry up to n copies of 0(2). 
In passing, it should be noted that if the random bond 
variables Jy for spin 1 and spin 2 are independent -which 
is the case Dotsenko and Dotsenko mainly studied,E3 crit- 
ical behavior is slightly different form the one below. 

Let us set g + = -<?_ = g M (§v = 0) and g A = 0. 
Then, the scaling equations for this specific model are 
given by 



dc/M 
dl 
dg 

n = ~ 99M - 



-9m+99m, 



(4.2) 



In the following subsections, we study the 5 = case and 
g 7^ case separately. 



A. Non-interacting case 

This subsection is devoted to the case g = 0, which cor- 
respon ds to the random bond Ising model. As a solution 
of Eq. (4.2), we have 



The scaling equations in the case, j / have been 
solved by Dotsenko and DotsenkorZl It may be conve- 
nient t o d erive first the relationship between c/m and g: 
From (O) 



9m(1) = 9(1) In I ~7j) ) > 



9 



where 



9* = 9o exp 



( 9M0 \ 


2ffo 


/ .9a/o\ 




= exp 


\2goJ 


V 9o J 


7T 



(4.6) 



(4.7) 



This constant is a fixed point value of g for negative go, 
as we shall see momentarily. These equations lead to 



(4.8) 



M l )J \9o 



where li(a;) is the logarithmic integral function. In what 
follows, we solve the equations separately for positive and 
negative cases. 



9m 



9mo 



1 + 9mo(1 - h)' 



(4.3) 



where cjmo is the initial coupl ing co nstan t at I = lg. Sub- 
stituting this into Eqs. (2.11) and ( 2.12| ), we have 



In - = I - l + i In [1 + g M0 (I - l )} 
~l-lo + -hx(!-l ), 



(4.4) 



and together with (2.13), we obtain the following DOS 
p(*0~fln(A). (4.5) 

The nonlinear sigma model for the class D has 
predictecSE2l |lnU| behavior of DOS. Actually, on in- 
termediate e nerg y scales, the | \nE\ factor dominates the 
DOS in Eq. Q. 



B. Interacting case 

If the interaction is finite, there are two fixed points 
depending on the sign of the interaction go as follows: 

1. <7Af(0 — * and g(l) for 50 > 0. 

2. ~g M {l) -> and g(l) -> g*{< 0) for g < 0. 

This indicates that the Ising fixed point or, in other 
words, the class D is stable against the positive inter- 
action. In the case of the negative interaction, however, 
it flows to fixed points depending on go, away from the 
class D. 



1. The go > case 

As mentioned above, g*/g(l) — > +00 in this case. By 
the use of the asymptotic expansion of the logarithmic 
integral function 

lia; ~ x (7— — I 5 h...|, for x -> +00, (4.9) 

Vina; In 1 / 

we can obtain the following asymptotic expression of the 
coupling constants 



g(l) ^ 9 * (1 - lnlnx + _L 
xlnx V \nx \nx 



x \ mx 



(4.10) 



with x = g*l. Eq. ( [2.12 ) is computed asymptotically as 
In^I + ilnl + lnlnl + of^jJ. (4.11) 



This equation, together with Eq. (2.13), leads to the 
following DOS, 



(4.12) 



As compared with the previous result (4.5), extra 
In I In E I factor shows up due to the flow effect of the 
coupling constant g. It should be noted, however, that 
this fixed point belongs to the class D. 



4 



2. The go <0 case 

The result obtained so far tells that universality class 
D is stable against the positive interaction. However, for 
the negative interaction, g(l) flows to nonzero g*(< 0) 
and a new universality class appears. Note that < 
g*/g(l) < 1, then we use the following formula for the 
logarithmic integral function, 



lix ~ 7 + In \x — 1| 



x - 1 



. , for x~l, (4.13) 



where 7 is Euler's constant. Now we have 

~9* 



9 ~ 
9m 



1 - e x 
9* 



1 



ln(l-e*), 



(4.14) 



where x is the same definition as the one in Eq. (4.10) 
Therefore, Eq. (2.12) is computed as 



In- 



A 

E 



zl + 0(e-\' 9 ^ 1 ), 



1,9* 1 1 1 1 ( 9mo 

1 + T = 1 ->I° P (-W 



and we reach 



P(E) 



(2-z)/z 



(4.15) 



(4.16) 



The power-law behavior of the DOS tells us that that the 
criticality is controlled by the interaction, and the fixed 
point does not belong to the class D. 



C. Spin correlation functions 

So far we have calculated the scaling dimension of the 
ensemble-averaged two-point correlation function of the 
fcrmion fields, and DOS as a result. We will calculate 
in this subsection spin correlation functions of the disor- 
dered Ashkin- Teller model. 

With respect to the Majorana fermions x defined as 



4 



1 



(XlL ~ «X2i) 



V2 

i>L = "y= (XlL + iX2L) ■ 



(4.17) 



and similar for R-moving fermion, the random mass term 
and the interaction term can be written as 



Hm = iM {xirXil + X1RX2L) , 

'Hint = 2gXlRXlLX2RX2L- 



(4.18) 



The former equation tells us that the energy opera- 
tor e(z,z), which couples with the mass, is e(z,z) — 



£i{z,z) +e 2 (z,z) with Ej(z,z) = iXjR(z)xjL(z)- The 
operator product expansion of the energy operator and 
the s pin field, 04 (z, z), which is a scaling field of ctji in 
Eq. (U, di3t3 



i XjRX 3 L(z,z)a jl (w,w) 



1 



4-7T \Z 



r(jj{w,w). (4.19) 



Note our normalization of the coordinate z in Eq. (2.S) 
followed by 2ir. 

To calculate {a 1 {L)a 1 (0)) N and {<ti<t 2 {L) a^iQ)) 1 ^, 
let us consider the following two kinds of operators 



0[ N \z,z)=Y[a la (z,z) 



a=l 
N 



7 12 



( Z : Z ) = II a la(z, z)(T 2a (z, z). 



(4.20) 



Q = l 



We assume that there are enough num ber o f replicas 
n > N and compute t he O PE between (4.20) and four 
fermi operators in Eq. (2.7) using the basic OPE between 
Majorana fermions and the spin field ( 4.19Q . It turns out 
that the results depend on TV only, and we obtain the fol- 
lowing anomalous dimensions of the fields 0\ and O12, 
respectively, 



7 w 

72° 



N N(N-1)_ 



16 



-9m, 



2N 2N{2N~1) N 
-9M ~ -5-9- 



16 



(4.21) 



It should be noted that j[ N ^ is completely the same 
as the one for t he ra ndom-bond Ising model.liS Using 
Eqs. (|4.irj|) and (4.14), we have the correlation function 



(<Ti(I/)<7i(0)) Ar at large distances, 



N(N-l)/8 




(.90 > 0) 
(50 < 0). 



(4.22) 



As expected from the analysis in the last subsection, the 
correlation function of this type is not affected by the pos- 
itive interaction. In the case of the negative one, however, 
the logarithmic correction do not remain and correlation 
function is just the one for the pure Ising model. Con- 
trary to 7^^ , since 7^2 depends on g explicitly, it follows 
that the interaction effects are more manifest. Actually, 
we obtain 



{a l a 2 {L)u 1 a 2 {Q)) N 



5 



'Lo 



N/2 



In- 



L 



]V(2JV-l)/4 



In In — 

L 



N/4 



JV(l-§*/2)/2 



(.90 > 0) 

(.90 < 0). 

(4.23) 



Namely, due to the interaction, extra lnlnL factor shows 
up in the former case, and the exponent depends mani- 
festly on the cou pling constant of the interaction through 
the relation ([O]) in the latter case. 



V. THE RANDOM SCALAR POTENTIAL 



Setting g + = g_ = gv as well as gA = we have 



dgv ~2 
-^ f =9v + 99v, 

dg 
dl 



99v- 



As discussed by Ludwig et.al.J3 the case with random 
scalar potential only (g = 0) belong to the class All 
because of time-reversal symmetry]!! It follows that the 
model is expected to flow to the fixed point of the sym- 
plcctic nonlinear sigma model, although in the present 
weak-coupling approach, gy is marginally relevant and 
flows to a strong coupling regime. With finite go > 0, 
Eq. (5.1) tells us that gv and g also flows to infinity. In 
the case go < 0, however, the model has a fixed point 
gv — > and g — > g*, where g* is specified momentarily. 

In this section, we examine only the negative go case, 
and follow a similar calculation to the previous section 
of the random mass model. First of all, we note the re- 
lationship 



~ 9v{l) = ~~ 9ln {W) 

9vo \ 
9o J 



g* = .9o exp 



(5.2) 



which correspond to Eqs. ( |4.6| ) and fl4.7|) . Therefore, the 
asymptotic solutions for large I are given by 



z = 1 



= 1 



ISOI 

7T 



exp 



9vo 
2|ffo| 



(5.5) 



It turns out that negative interaction yields nontrivial 
fixed point dependent on the strength of the interaction, 
which does not belong to AIL 



VI. THE RANDOM POTENTIALS WITH 
CHIRAL SYMMETRY 

So far we ha ve studied some fixed points of the RG 
equations ( 2.10 ). If we allow non-hermitian Hamiltonian, 
there is-still another fixed point with <?+ = 0. Guruswamy 
et. al.EI have proposed a model with a random imagi- 
nary vector potential, a random imaginary scalar poten- 
tial and a random real mass term, as a fermion model 
equivalent to the random phase sine-Gordon modelliSEj 
In this section we study how interaction affects the crit- 
ical behavior. Let us add the interaction (2.3) to the 
Hamiltonian which Guruswamy et. al. have studied 



( 51 ) H = ipi^n (d^ + Ap) ip + Mtpa3ip + iVipip + "Hint- (6.1) 



This is basically the same model as the random phase 
sine-Gordon model with arbitrary boson coupling K.E3 
The RG equations of this model arc given by setting 
9A -* ~9A, g+ = (gv = -9m < 0), and = -2g M -> 

-.9-, 



dgA 
dl 
dg~ 
~dT 
dg_ = 
dl 



gg- 

-99- 



(6.2) 



This model may not be suitable for the computation 
of the DOS because of its non-hermiticity. A simple ex- 
tension, however, to a two Dirac fermion model yields a 
hermitian Hamiltonian with the same RG equations in 
the following way; 



H12 = T~ti + Ti-2, 



(6.3) 



where Tij is the Hamiltonian ( p.l[ ) for the fermion ipj and 
ipj. Next, make the transformation, 



9 = 

gv 



1 + es* 1 ' 
9* 



■]n(l + e s ' 1 ). 



1 + eS' 1 

Using these, we finally end up with 



P{E) 



A 



(2-z)/z 



(5.3) 



(5.4) 



t_2R ->■ ^2fl, 1p2R -> Tp 2 Ri 
'4>2L i>lL, ^2L IplL, 



(6.4) 



which yields no change to the partition function up to 
some constants, and, at the same time, the transforma- 
tion as well 



where 



Ai -> A 2 , A 2 -> -Ax, 



(6.5) 



G 



which keeps the probability distribution (2.6) invariant. 
Then, these transformation yields 

H12 -» ^i«7ju (dp - iA^) 4>i + ^2*7^ (df* + iA n) V>2 
+Bip 1 a 3 ^ 2 + B1P2V31P1 
+2g(Ji R J2 L + J 2 RJi L ), (6.6) 



A. Non-interacting case 



To begin p4th, we review the results known so far for 
<7o = caseJla'S The scaling equations for c/a and <?_ in 
( |6.2| ) are easy to integrate, giving the flow of the coupling 
constants, 



where B is defined here as B — M + iV and JjR(L) — 
^jRiL^jRiL) denotes the current of the jth fermion. 
This Hamiltonian is actually hermitian, and probabil- 
ity distribution of is the same as in ( |2.6| ), and 

B(x)B(y)=g_5(x-y). 

Without the interaction (g = 0) in Eq. ( |6.6| ), this 
Hamiltonian is just the continuum limit of the random 
hopping fermion model on a square lattice with 7r-flux per 
plaquettej3 or the one on a honeycomb lattice. It should 
be noted that the interaction term in the Dirac Hamil- 
tonian ( |6.6| ) is not due to interactions, but due to e.g^. 
annealed randomness, of the original lattice models.E2l 
Actually, annealed randomness for quenched random sys- 
tems have turned out-tp play an important role in metal- 
insulator transitionsE3 With n replicas, this extended 
model belongs to the class BDI when g — 0, and has a 
global U(2n) symmetry, which include chiral U(n)xU(n) 
symmetry. To see this, it may be convenient to write the 
above Hamiltonian with g = as follows; 

W12 = V R 2(d- z + A^r + ~* L 2(d z + A Z )^ L 

-iB^ L ^R + iB^R^ L , (6.7) 

where ^ R = {ipi R ,ip 2R )> *k = (^lfl; ^2r), and similar 
for and ^l- It is readily seen that the replicated 
model is invariant under the transformation 



MO = 3-o, MO = 9A0 + g-o(l - la) 



Vr 



^R 



>U*R, 



(6.8) 



with 2n x 2n unitary matrix U. Of course, the inter- 
action breaks this symmetry, as in the previous cases. 
Interestingly, this Hamiltonian has another symmetry. 



*7V. 



with 



B -> Be 



2 1 1) 



B 



e ie V R , 



'B. 



(6-9) 



(6.10) 



This symmetry as well as ( |6.8D allows us to cetermine 
the gu-dependence of the correlation functions.EScI 

We study this model perturbed by — z^^! +i/j 2 i/> 2 ) to 
derive the scaling function of the DOS. The anomalous 
dimension of this operator is 



7 



1 



MO , MO 



4 ' 2 ' 

and the scaling equations are the same as Eqs. (|6.2| 



(6.11) 



(6.12) 



where §ao and g~o are the initial coupl ing constants. 
Substituting theses solutions to Eq. ( 2.12| ), we obtain 



(6.13) 



where 



n 2 



9ao g-v_ _ 1 
4 2 

9mo 



gAo + gno 



7T 



(6.14) 



Combining this equation with Eq. ( [2.13] ), we reach the 
famous expression of the DOS -for the class BDI, which 
Gade derived for the first timer! 



/ x A 

P(E) ~ exp 




(6.15) 



This formula shows that the DOS diverges towards the 
zero energy. However, the crossover energy scale under 
which the singularity of the DOS emerges is estimated as 



E r 



Ke~ z2/{2v) ~ A exp 



2.9mo 



(6.16) 



or in terms of the length scale, the crossover length is 
given by L cr ~ L exp (27r 2 /g^ ) , followed from Eq. 
Q6.12 ). Namely, the crossover energy (length) is expo- 
nentially small (large) for weak randomness. Considering 
this, we expect 



P(E) 



A 

E 



exp 




In-) (E«E cr ) 

y e 



f) 



(2-z)/z 



(6.17) 



(E » E cr ). 



Previous numerical studiesB observed a power-law behav- 
ior rather than the singular divergence even quite near 
the zero energy. In the regime of their study, E cr may be 
too small to yield the divergent DOS, or the system size 
is much smaller than the crossover length. Actually, Ryu 
and Hatsugai have recently been successful in observing 
the divergent DOS by studying quite large systems.EJ 
The DO S cal culated by them may be described directly 
byEq. floTSl). 
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B. Interacting case 



2. The case c = Q 



Next consider the case with finite interaction go =/= 0. 
First, it should be noted that g_ + g is conserved under 
the scale transformation, 

9- + 9 = 9-0 + 9o = c (6.18) 

Therefore, only two equations among three in Eqs. ( |6.2| ) 
are essentially independent. By using c, the equation for 
<?_ can be converted into 



dg- 
dl 



-g-(g- - c). 

It is easy to solve this equation, 

MO 



1 + (g /g-o)e-< l - l o) ' 
for nonzero c, whereas for c = 0, 

9-0 



M0 = 



1 + 9-oQ - l o)' 



(6.19) 



(6.20) 



(6.21) 



Directly from these solutions or by analyzing the beta 
functions in Eq. ( |6.2j ), it is not difficult to observe that 
there are three kinds of fixed points. 

1. g~(l) — > c (= g_ + g ) and g(l) — > for c > 0. 

2. g_(l) -> and <?(/) -> for c = 0. 

3. <7_(Z) -> and g(l) -> c for c < 0. 

In the following subsections, we study the cases sepa- 
rately. 

1. The case c > 

As g(l) — > 0, this fixed point is just the same as that 
of the non-interac ting case. The flow of <7 a(0 is easy to 
integrate in Eq. (6.2) by the use of Eq. (6.20), and for 
large I the coupling constants are approximately given by 



M0 ~ c ' 

M0 ~ c 2 (l - l ) + gAo - cln 



9o 



Substituting these into Eq. (2.12), we have 



]n±~z'(l-l ) + ^(l-l ) 2 , 



g . (6.22) 



(6.23) 



with 



z + 



9o ~ {9 mo + 9a) In 



( 9o 



\9mo 



c 2 /4 



go + 9 mo 



(6.24) 



where z is d efined in Eq. (6.14). Above equation leads 
to Eq. ( 6.1 5| ) with modified parameters y' and z' . 



We expect, at first sight, that this case also flows to 
the same fixed point as above, since g(l) — > 0. However, 
the integration of c/a gives 



M0 



9-0 



l+g-o(l-lo) 



+ 9-o +9 AO- 



(6.25) 



I n the previous cases, the I term of §a in Eqs. ( |6.12j ) and 
( |6.22 ) is responsible for the divergent DOS at the band 
center, which is missing in the present case. Now we have 



In A 

E 



z(l-l ) + -ln(l-lo), 



with 



1 



gAo + 9-0 



2.9 AO + 9 mo 
2^ ' 



This leads to 

P {E)r 



(2-z)/z 



In 



X/2z 



(6.26) 



(6.27) 



(6.28) 



The DOS in this case shows a power-law behavior with a 
logarithmic correction. Although this fixed point should 
be in the class BDI, the behavior of the DOS is quite 
different. 



3. The case c < 

As we have seen in the last two subsections, the class 
BDI is stable against the interaction in the case c > 0, 
although DOS is different from that of BDI for c = 0. 
Contrary to this, the present case leads to a different 
fixed point which does not belong to BDI any longer. 
Asymptotically for large I the flows of the coupling con- 
stants are computed as 



M0 
M0 



0, 

<?A0 



go - c In 



.90 

9-o 



0{e 



-|c|h 



Therefore we obtain 



4 



z(l - l ), 



(6.29) 



(6.30) 



where 



1 



9 AO 



9o 1 , 
-cm 



4 4 

2g Ao + gMo 



1 



2vr 



+ 



.9-0 

\9mo+9o\ 
2tt 



In 



9o 



9mo 



This equation leads to 
p{E)~ 



(2-z)/z 



(6.31) 



(6.32) 
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VII. SUMMARY AND DISCUSSIONS 



In this paper, we have investigated the interaction ef- 
fects on the disordered Dirac fermion model separately 
with a random vector potential, a random mass, a ran- 
dom scalar potential, and their special combination with 
chiral symmetry. They belong to the classes AIII, D, All 
and BDI, respectively, when the interaction is zero. It 
turns out that some fixed points are stable against the 
interaction in the sense that the coupling constant of the 
interaction flows to zero in the RG, whereas others are 
unstable and flow to fixed points which do not belong to 
those classes listed-above. Thus, it turns out that interac- 
tions play an important role in the symmetry properties 
of localization problems. 

So far we have studied several specific fixed points. In 
what follows, we would like to discuss briefly the case 
where all kinds of disorder are simultaneously present. 
As was shown by Ludwig et. al.,L3 the fixed line of the 
random vector potential as well as the fixed point of a 
random mass are unstable and flow to unknown strong- 
coupling regime on generic initial conditions. In the 
present case with the positive interaction, the situation 
is quite similar to the case of non-interactin g cas e. Nu- 
merical calculation of the scaling equations ( |2.10| ) shows 
that the coupling constants flow to infinity. The neg- 
ative interaction, however, can give rise to a different 
situation. Below a certain critical g = g c < 0, which 
may depend on the other disorder strength, the solution 
of the scaling equations flow to a fixed point, which is 
quite similar to the fixed line of the model with the ran- 
dom vector potential and the interaction only, discussed 
in Sec. III. Namely, after the flow, qa and g become 
finite constants whereas g± — > 0. Therefore, we could 
claim that the fixed line formed by the random vector 
potential and the marginal interaction is one of generic 
critical points for disordered Dirac fermion model with 
randomness and interaction. Other fixed points are in a 
strong coupling regime and are not accessible by a weak- 
coupling approach. 

As it has turned out that an interaction plays a role in a 
simple disordered model like Dirac fermions examined in 
this paper, it is quite interesting to study an interplay be- 
tween interactions and randomness for more complicated 
systems, especially those belonging to the new universal- 
ity classes. 
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